NAVIGATION PROCESSING FOR A SATELLITE POSITIONING SYSTEM RECEIVER 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application claims priority under 35 U. S. C. § 119(e) of United States Provisional 
5 Patent Application No. 60/228,965, filed August 29, 2000, entided "NAVIGATION 

PROCESSING FOR A SATELLITE POSITIONING SYSTEM RECEIVER," by Keith J. Brodie 
et al., which application is incorporated by reference herein. 



10 BACKGROUND OF THE INVENTION 

1. Field of the Invention. 

The present invention relates in general to receiver operation in satellite positioning systems 
(SATPS) such as the Global Positioning System (GPS), GLONASS, and the proposed Galileo 
system, and in particular to navigation processing for a SATPS receiver. 

15 

2. Description of the Related Art. 

Navigation processing for SATPS receivers usually comprises processing pseudoranges and 
Doppler frequency shift measurements to determine position and velocity of the receiver. However, 
this approach is sometimes unable to distinguish between alternate (valid and invalid) three-satellite 
20 solutions. Examples of such processing techniques are found in U.S. Patent Nos. 5,416,712 by 
Geier, US 5,883,595 by Colley, and US 5,590,043 by McBurney, all of which are incorporated by 
reference herein. 
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SUMMARY OF THE INVENTION 

t 

To minimize the limitations in the prior art, and to minimize other limitations that will 
become apparent upon reading and understanding the present specification, the present invention 
discloses a method and apparatus for navigation processing in a satellite positioning system receiver, 
including methods and apparatus for generating a least-squares fix, checking the validity of 
measurements, and processing measurements in a low-power mode which involves repeated on/off 
cycles of the reference oscillator used to make the measurements. 

A method in accordance with the present invention comprises separating the three SATPS 
satellites into a first pair and a second pair, constructing a primary solution and an alternate solution, 
wherein the primary solution and the alternate solution satisfy the measurement constraints, 
computing a Doppler difference estimate for the primary solution and a Doppler difference estimate 
for the alternate solution, computing Doppler difference residuals for the first pair and the second 
pair of SATPS satellites, and comparing the Doppler difference residuals for said primary and 
alternate solutions to determine a valid solution. Typically, the computing of a Doppler difference 
residuals comprises differencing a measured Doppler difference from an estimated Doppler 
difference for the first pair and the second pair of SATPS satellites. Typical determination of a valid 
solution comprises comparing the difference between the Doppler difference residuals to a 
predetermined number. Usually when the Doppler difference residuals exceed the predetermined 
number, the alternate solution is selected as the valid solution. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

t 

Referring now to the drawings in which like reference numbers represent corresponding 
parts throughout: 

FIG. 1 illustrates the geometry of alternate solutions for a three satellite fix; 

FIG. 2 illustrates the geometry of the projection vector used to compute the initial point for 
the least-squares iteration for the alternate fix solution; 

FIG. 3 is the flow chart for the least-squares fix decisions on posting computed fixes and 
testing for alternate soltutions; 

FIG. 4A and 4B are the flow chart for initial measurement validity checking with an adaptive 
form of a signal power level test. 

FIG 5. is the flow chart for the decision to use, or not use DGPS corrections to GPS 
measurements. 

FIG 6A,B, and C are the degraded mode logic for SATPS receiver 
FIG 7 is the Receiver Manager to Navigation Library Interface 

FIG 8A through 8G is the flowchart for the main one-hertz navigation processing loop. 
FIG 9. Shows the navigation mode decision with an overview of the processing in each 

mode. 

FIG 10. is a flowchart of the mode decision logic for the main navigation loop. 



siRF.tnusui 



DETAILED DESCRIPTION OF THE DRAWINGS 

t 

In the following description of the preferred embodiment, reference is made to the 
accompanying drawings which form a part hereof, and in which is shown by way of illustration a 
specific embodiment in which the invention may be practiced. It is to be understood that other 
embodiments may be utilized and structural changes may be made without departing from the scope 
of the present invention. 

Least-Squares Navigation Fix 

A least-squares iterative solution for position and clock bias with four or more satellites is well 
known. A description of the basic least-squares algorithm for GPS pseudoranges can be found in 
Understanding GPS: Principles and Applications by E. Kaplan (ed.), Artech House, 1996, pp. 43-47 
which is herein incorporated by reference. A three-satellite solution can also be computed by using 
altitude as a fourth measurement. If the altitude is known it may be used, if not, a default altitude 
may be assumed for receivers designed for terrestrial use. Altitude may be known from user input 
or from a saved altitude from a prior period of navigation with more than three satellites. 
There are usually two solutions to the three-satellite problem. A diagram of the situation is shown 
in Figure 1 . The sphere is a representation of a constant altitude surface over, or under, the Earth's 
surface. In other words, every point on the sphere satisfies the altitude constraint. The Earth's 
surface is not spherical of course, but close enough so that the diagram illustrates one instance of 
alternate three satellite solutions. The two points identified on the sphere satisfy the pseudorange 
constraints with different values for user clock bias, or we could say, equivalently, that the difference 
in the ranges to the satellites is the same for the two alternative fixes. 
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In order to safely report and use a three satellite least-squares solution, we must be able to 
determine if the solution we have computed is the correct one from the two possibilities. We have 
additional data to assist in checking the validity of the solution. The difference in doppler 
frequencies of the satellites in track observed, versus that predicted from the known satellite orbits 
and the computed solution is a test on the validity of a computed root. The doppler difference 
residuals are used to remove the error in the receiver's clock frequency from the test, so the test has 
finer resolution and is generally applicable under a wider range of circumstances than a test on 
absolute doppler. 

The doppler difference residual is computed by differencing the observed doppler 
frequencies between two satellites in track. This difference is independent of the frequency error in 
the SATPS receiver's oscillator by which the received signal's doppler is measure. An estimate of 
this same doppler difference is made by assuming the SATPS receiver is at the location of a 
proposed fix, and computing the doppler for the known satellite orbits by computing the dot- 
product of their earth-fixed frame velocity vectors with the line-of-sight vectors from the fix 
location to the satellite positions. The doppler difference residual is then the measured difference 
between doppler frequencies for two satellites minus the estimated difference between dopier 
frequencies for two satellites. The sum of the absolute value of the doppler differences is then 
checked against thresholds for reporting and testing alternate roots. The sum of the squared 
residuals could be used equally well with appropriate changes in the threshold values used to make 
the reporting and testing decisions. 

Unfortunately, the differenced doppler residuals between two alternate roots are not always 
large enough to guarantee proper root selection. Any velocity of the receiver, in addition, is another 
unknown in the doppler difference residuals which will limit the ability to decide between roots. In 
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order to safely use one of the roots, the second root should be computed and checked. If intial user 
velocity is known, or sensor data such as a wheel speed sensor and a compass or heading estimate 
mantained by a gyro is available, the speed data can be included in the relative velocity computation 
used for generating the estimates. 
5 A fast method to compute the second root has been found using a projection of the first 

root through the plane of the satellites. This method results in an approximate location for the 
second root, which is used to initialize the least-squares iterator so that it will close on the alternate 
root. The position of the three satellites used in the solution is used to form a plane. A line is 
drawn normal to this plane and intersecting the first root found on the surface of the constant 
10 altitude sphere. This geometric construction is shown in Figure 2. There are cases where the plane 
of the satellites does not intersect the sphere, nevertheless there may be two points on the normal to 
the plane that intersect the sphere. The projection point is computed from the initial root, the 
position of the satellites, and the radius of the constant altitude sphere using the following method. 
S x , S 2 , S3 : Satellite position vectors 

15 R x : First root computed from least-squares iterator 
Satellite position vector differences D x = S x - 5 2 
D 2 =S 2 -S, 

Compute a projection vector, normal to the plane defined by the three satellites. 

P = D { xD 2 

20 Compute the unit-normal projection vector. 




The form of the second, projected root is assumed to be 
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R 2 =R { +aP 

r 

where a is an unknown constant to be solved for. Since the second root is also required to lie on 
the constant altitude sphere: 



\R 2 \ = \R l \ = \R l +a-P\ 

This results in a quadratic for a with one non-trivial root: 

a = -2R } P 

The projected root is then 

R 2 =R^2(R r P)P 

R 2 is then used as the initial position estimate for the least-squares estimator, which will normally 
close on the alternate root. Occasionally this method may fail to find an alternate root to check. In 
this case, the three satellite fix is rejected, rather than reporting a fix which may be substantially in 
error. In some cases, the doppler difference residuals observed at the first root and the alternate 
root may not differ by a substantial amount. In this case, again lacking clear information to choose, 
the fix is rejected rather than reporting a fix which may be substantially in error. 

Measurement Validity Testing 

There are many potential sources of error in SATPS system measurements. These include noise, 
tropospheric and ionospheric delays, multipath, errors in syncing to the CDMA code used to 
modulate the satellite transmissions, auto-correlations with the wrong portion of the code, cross- 
correlations with a different code from another satellite, interfering signals from the local 
environment or other satellite systems, frequency alias in the tracking loops, and others. 
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Errors such as receiver measurement noise, atmospheric correction errors, errors in the predicted 
satellite positions are well understood and modeled in the navigation processing. The other errors, 
such as code-sync errors, multipath, and auto- and cross-correlations are not well modeled, can be 
very large, and would severely affect navigation processing if measurements affected by these errors 
are not identified and removed from the set of measurements to be processed. Accordingly, a 
strong set of validity tests has been developed to insure the integrity of the measurements prior to 
processing. Those tests can largely be divided into those that depend on reasonable knowledge of 
the current navigation state, and those that are independent of the state estimate and that can 
therefore be applied in all conditions, whether prior navigation data exists or not. 

One of the novel test concepts utilized was the creation of an adaptive power mask. The 
basic power mask test is a check on the carrier to noise density, or C/N 0 that the signal was tracked 
at to insure a reasonable prospect that the measurement is valid. This basic test has been modified 
to make it adaptive in two ways, as shown in Fig 4. First, the threshold power level is adapted to the 
current navigation mode. If the SATPS receiver is navigating, there are a number of subsequent 
measurement validity checks that will be performed, so it is reasonable to lower the power threshold 
in the validity test to attempt to maintain the largest possible number of available good 
measurements. If, however, the receiver is not navigating, the power threshold is kept higher to 
prevent measurements with errors from getting into the initial solution. The second adaptation 
mechanism is to vary the form of the test between an average and minimum power test. The 
average power test is less restrictive, will allow more measurements, and has increased reliance on 
subsequent testing to find measurement errors. The minimum power test is more restrictive, 
requiring that all of the 100 millisecond bins used to generate the measurement have power levels 
above the threshold. Since the 100 millisecond estimate of C/N„ is noisy, this test will result in the 
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occasional rejection of measurements from signals near the threshold from C/N„ estimate noise 

f 

alone, but this is a reasonable trade-off to make against the risk of a measurement with large errors 
being used in the initial navigation fix. 

Following the loop defined in Fig. 4 the set of valid pseudoranges above the power 
threshold has been identified. The set of valid delta-pseudoranges are identified as those coming 
from the same measurement channels as the valid pseudoranges with the additional criteria that the 
coherent integration time for the carrier signal be non-zero. 

A test is performed on the absolute value of the differences in the delta-pseudoranges. This 
test is independent of the current navigation mode and can be applied without any knowledge of the 
SATPS receiver's position, the time, or the oscillator frequency error. In this test, the difference 
between to delta-pseudorange measurements is compared to the absolute physically possible limit 
for a terrestrial user given the constellation of satellites in use by the SATPS. In the current GPS 
orbits, for example, the maximum doppler to a stationary terrestrial user is about 800 m/s, so that 
the maximum doppler difference, observed between two satellites tracked simultaneously, must be 
less than approximately 1600 m/s. User velocity, taken to be a small number relative to these 
satellite dopplers, does not change the limit substantially because a significant user velocity towards 
one satellite would be substantially away from another satellite which was far enough away (in 
azimuth to the receiver) to be close to the limit. Thus all available differences of valid delta- 
pseudoranges are computed and compared to the 1600 m/s threshold. For any pair that violates the 
threshold, a counter for each of the participants in the offending difference is advanced. At the end 
of the check, a single bad measurement may have caused large differences in multiple pairs, but the 
counters will identify the culprit assuming sufficient satellites are available. If any measurement is 
involved in more than one bad difference, it is rejected, and if no bad measurement, for example, 
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with two valid delta-pseud 9 ranges with a bad difference, they are both rejected. The fault has been 
identified but not isolated. 

A completely analogous check is made on the pseudorange measurements. For GPS, the 
physical range from a satellite to a SATPS which can receive the signal is limited to a range from 
about 20 million meters for a satellite overhead to about 26 million meters on the horizon. Without 
reference to the navigation state, the time of day, or any other state estimate, an absolute check < 
pseudorange measurements can be performed. No two GPS pseudoranges recorded at the same 
time for a terrestrial STAPS receiver can differ by more than 6 million meters. Counters for bad 
differences are kept as in the delta-pseudorange difference check to identify faulty measurements 
from in a set of difference checks. 
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Low-Po wer Measurement Processing 

In low-power processing the reference oscillator is turned-off between measurement update: 
to save power. As a result, substantial oscillator dynamics, due to heating, are present in each set of 
measurements. These reference oscillator dynamics adversely affect normal processing of the delta- 
pseudorange measurement in a Kalman filter because of their unmodeled character. The filter must 
assign these unmodeled accelerations to either the clock drift term or velocity terms in it's state 
update, and errors are inevitable. In order to work around this problem, differenced delta- 
pseudorange measurements are processed instead of the whole delta-pseudorange. This limits the 
observability of clock drift in the filter and results in substantially better velocity estimates than 
would be otherwise obtainable. 
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Conclusion 

A novel method for navigation processing in a satellite positioning system receiver has been 
disclosed. 

A method in accordance with the present invention comprises separating the three SATPS 
satellites into a first pair and a second pair, constructing a primary solution and an alternate solution, 
wherein the primary solution and the alternate solution satisfy the measurement constraints, 
computing a Doppler difference estimate for the primary solution and a Doppler difference estimate 
for the alternate solution, computing Doppler difference residuals for the first pair and the second 
pair of SATPS satellites, and comparing the Doppler difference residuals for said primary and 
alternate solutions to determine a valid solution. Typically, the computing of a Doppler difference 
residuals comprises differencing a measured Doppler difference from an estimated Doppler 
difference for the first pair and the second pair of SATPS satellites. Typical determination of a valid 
solution comprises comparing the difference between the Doppler difference residuals to a 
predetermined number. Usually when the Doppler difference residuals exceed the predetermined 
number, the alternate solution is selected as the valid solution. 

The foregoing description of the preferred embodiment of the invention has been presented 
for the purposes of illustration and description. It is not intended to be exhaustive or to limit the 
invention to the precise form disclosed. Many modifications and variations are possible in light of the 
above teaching. It is intended that the scope of the invention not be limited by this detailed description, 
but by the claims appended hereto. 



SiRF.mUSUl 



